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A test of the bosonic spinon theory for the triangular antiferro- 
magnet spectrum 
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PACS 75. 10. Jm - Quantized spin models, including quantum spin frustration 

Abstract. - We compute the dynamical structure factor of the spin-i triangular Heisenberg 
model using the mean field Schwinger boson theory. We find that a reconstructed dispersion, 
resulting from a non trivial redistribution of the spectral weight, agrees quite well with the spin 
excitation spectrum recently found with series expansions. In particular, we recover the strong 
renormalization with respect to linear spin wave theory along with the appearance of roton-like 
minima. Furthermore, near the roton-like minima the contribution of the two spinon continuum 
to the static structure factor is about 40% of the total weight. By computing the density-density 
dynamical structure factor, we identify an unphysical weak signal of the spin excitation spectrum 
with the relaxation of the local constraint of the Schwinger bosons at the mean field level. Based 
on the accurate description obtained for the static and dynamic ground state properties, we argue 
that the bosonic spinon theory should be considered seriously as a valid alternative to interpret 
the physics of the triangular Heisenberg model. 
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Introduction. — During a long time the magnetic 
ground state of the spin-i triangular Heisenberg model 
(THM) has attracted the attention of many researchers, 
due to the possible realization of the resonating valence 
bond (RVB) ground state proposed by P. W. Ander- 
son in 1973 [I]. The revival of the RVB theory for the 
cuprates [2] prompted the investigations of quantum dis- 
ordered ground states within large N theories where the 
Heisenberg interaction is naturally written in terms of 
singlet bond operators and fractional spin-i excitations 
with bosonic or fermionic character [3J. The fermionic 
i version leads to exotic disordered ground states H] while 
the bosonic one allows to describe disordered and ordered 
ground states [5] by relating the magnetization with the 
condensation of bosons [6]. For this case, using gauge 
field theoretical arguments, it has been conjectured that, 
when short range spiral correlations are present in the dis- 
ordered phases, the bosonic spinons would be in a decon- 
fined regime [5]. Therefore, a broad two spinon continuum 
is expected in the spin excitation spectrum. 

From the numerical side, instead, thanks to the enor- 
mous effort of the community to develop unbiased tech- 
niques [7HTU]. it has been firmly established that the 
ground state of the spin-i THM is a robust 120° Neel 
order. These numerical results precluded the fermionic 



version of the RVB theory, giving support to both the lin- 
ear spin wave theory (LSWT) and the bosonic version of 
the RVB theory, namely the Schwinger boson theory. In 
fact, both theories agree quite well with numerical results 
on finite size systems [TT1IT2"] , although for spiral phases the 
singlet structure of the mean field Schwinger bosons theory 
does not recover the spin wave dispersion relation in the 
large s limit [13]. Consequently, linear spin wave theory 
seemed to capture the quantum and semiclassical features 
expected for a 120° Neel ground state of the THM. How- 
ever, recent series expansions studies [TUJ[T3] challenged 
LSWT, showing that for s = \ the functional form of the 
dispersion relation differs considerably (points of fig. [2]) 
from that of LSWT (solid line of fig. [2]). In particular, 
it was observed a strong downward renormalization of the 
high energy part of the spectrum along with the appear- 
ance of roton-like minima at the midpoints of edges of the 
hexagonal Brillouin zone (BZ) (B and D points of the inset 
of fig. [I]). The authors argued that the differences with 
LSWT could be attributed, probably, to the presence of 
fermionic spinon excitations. Nevertheless, further spin 
wave studies [H] showed that, to first order in 1/s, there 
appear non trivial corrections to the linear spin wave dis- 
persion due to the non collinearity of the ground state, 
giving a fairly accurate description of the series expansion 
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results. However, magnons are not well defined for an am- 
ple region of the BZ [TB]. Another question, regarding the 
spectrum of the THM, is the nature of the multiparticle 
continuum above the one magnon states. For instance, it is 
believed that the broad multiparticle continuum measured 
in the CsiCuCl^ compound is better described by an in- 
teracting spinon picture than a magnon one [171 . In this 
sense, given accurate predictions of the Schwinger boson 
theory for the static ground state properties of the THM 
[11] , it is important to investigate whether the anomalous 
features of the spectrum found with series expansions can 
be captured, or not, by this alternative theory that natu- 
rally incorporates fractional spin-^ excitations. 

In the present article we investigate the validity of the 
bosonic spinon theory to interpret the spin excitation spec- 
trum of the spin-i THM. Our main finding is that the 
mean field Schwinger boson theory (intensity curves of fig. 
[2]), based on the two singlet operator scheme [T5], repro- 
duces qualitatively and quantitatively quite well the re- 
cent series expansions results. By computing the dynami- 
cal structure factor, we remarkably find that the expected 
spin excitation spectrum is recovered by a reconstruction 
resulting from a non trivial redistribution of the spectral 
weight located at the spinonic branches shifted by ±§, 
where Q = (|tt, 0) is the magnetic wave vector. By com- 
puting the density-density dynamical structure factor, we 
were able to identify, at the mean field level, the remnant 
weaker signal of the spectrum with the relaxation of the 
local constraint of the number of bosons. We also discuss 
the validity of the alternative mean field decoupling based 
on one singlet operator scheme. 

Mean field Schwinger bosons approximation. 

In the Schwinger boson representation [3] the spin op- 
erators are expressed as S,; = ^bja b^, with the spinor 
b} = (SL; St. ) composed by the bosonic operators St 



where 



and St,, and a— (a x , a y , <r z ) the Pauli matrices. To ful- 
fil the spin algebra the constraint of 2s bosons per site, 
So- ^leAc = ^s, must be imposed. Then, the spin-spin 
interaction of the Heisenberg Hamiltonian can be written 

as 



Sj-Sj — : BjjBij 



(1) 



where :: means normal order and the singlet bond op- 
erators are defined as = h a Ha^is an d B\- 



\ J2a b\J>j<j. We will briefly describe the main steps of 
the mean field while the details of the calculation can be 



found in our previous works |ll|H8j. Introducing a La- 
grange multiplier A to impose the local constraint on av- 
erage and performing a mean field decoupling of eq. (JlJ , 
such as Aij = (Aij) = (A^) and B tJ = {R tJ } = (Bt), the 
diagonalized mean field Hamiltonian results 



H 



MF 



= E t 



k 



[4 



-k4 a -H 



^gs = J$> k + AiV( S + -) 



is the ground state energy and 

fcfct = ^ = ^k = [(7k + A) 2 - (7k ) 2 ]*, 

is the spinon dispersion relation with geometrical factors, 
7]f = ^J^^iJacosk.J and y£ = ^J^2 S Ag sink. S, and 
with the sums going over all the vectors 6 connecting the 
first neighbours of a triangular lattice. The mean field 
parameters has been chosen real and satisfy the relations 
Bg = B-s and Ag = —A—g. The ground state wave function 
of Hmf can be written in a Jastrow form [6], 



|gs) = exp 



|0)&, 



(2) 



where |0}h represents the vacuum of Schwinger bosons 
and the odd pairing function is defined as = 
(/V)Sk/k elk(ri ~ r3) ' with ik = -v k /w k , and Bogoliubov 
coefficients u k = 7k ^ A )] ^ and v^ — i s<?n(7 k )[|(— 1 + 

7 * +A )]l. The singlet bond structure of eq. © guaran- 
tees the singlet behavior of |gs). Even if the Lieb-Mattis 
theorem cannot be applied to non bipartite lattices, the 
singlet character of the ground state for cluster sizes with 
an even number of sites N has been confirmed numerically 
7,8 . It should be noted, however, that |gs) is not a true 
RVB state because the constraint is only satisfied on av- 
erage. Furthermore, by solving the self consistent mean 
field equations at zero temperature, 



(3) 



A s = 


2iV^w k 


B s = 


\ UJi -cosk.(5 

2N ^ w k 

k K 


1 

S+ 2 = 


2N 4" w k ' 



it is found that as the system size N increases the singlet 
ground state |gs) develops 120° Neel correlations signalled 
by the minimum gap of the spinon dispersion located at 
±y, where Q = (|7T, 0) is the magnetic wave vector [T§] . 
As the spinon gap behaves as u , q ~ 1/N, for large system 
sizes the singular modes of eq. (J3J can be treated apart, 
analogously to a Bose condensation phenomena [6.. In 
particular, the local magnetization m(Q) can be derived 
from the last line of eq. ([3]) , yielding the relation [20] 



(7§ 



A) 2 



27V 



= S(Q) = ym 2 (Q), 



where S'(k)=^ R e lk R (gs|5o-S , R|gs) is the static structure 
factor. Formally, it can be shown that in the thermody- 
namic limit |gs) is degenerated with a manifold of Bose 
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Table 1: Energy and magnetization of the 120° Neel ground 
state of the spin-| Heisenberg antiferromagnet on the triangu- 
lar lattice as obtained with mean field Schwinger bosons within 
one [22] (A) and two [21] (AB) singlet scheme; Gaussian fluc- 
tuations [11] above the AB mean field (AB + Fluct), Quantum 
Monte Carlo [8] (QMC), linear spin wave theory (LSWT) and 
non linear spin wave theory (LSWT+l/s) [23] 





E/JN 


m 


A 


-0.7119 


0.328 


AB 


-0.5697 


0.275 


AB+Fluct 


-0.5533 




QMC 


-0.5458(1) 


0.205(1) 


LSWT 


-0.5388 


0.2387 


LSWT+1/.s 


-0.5434 


0.2497 



condensate ground states, each one corresponding to all 
the possible orientations, in spin space, of the 120° Neel 
order. In the Schwinger boson language the condensate 
of the up/down bosons at ±y and the normal fluids of 
bosons corresponds to the spiralling magnetization m(Q) 
and the zero point quantum fluctuations, respectively [6]. 
For the triangular lattice the present mean field approxi- 
mation [2T] gives a local magnetization m = 0.275. 

An alternative procedure, is to use the operator identity 
■.BjjBij-.+Alj Aij—S 2 , and write the spin-spin interaction 

|T]) in terms of the singlet operator A\j [3][221[M] : 

M i = -2ikA+ 5a - (4) 

Even if eqs. ([IJ and (@| are equivalent, the latter leads to a 
different mean field decoupling with parameters As and A 
[22 24]. In table[T]it is shown the values of the ground state 
energy and magnetization for the THM obtained with the 
two mean field Schwinger boson decouplings along with 
Gaussian fluctuations |11) . linear spin wave theory, non 
linear spin wave theory (LSWT+l/s) [53]; and quantum 
Monte Carlo [8 (QMC) results [25] ■ Even though it has 
not yet been calculated, we expect that Gaussian fluctua- 
tions above the mean field will reduce the magnetization, 
as has already been found for the spin stiffness in the THM 
[TT] . From table Q] it is seen that the two singlet scheme 
describe quantitatively better the static properties of the 
THM. 

Dynamical structure factor. 

Spin-spin correlation functions. We study the spec- 
trum through the dynamical structure factor at T = 0, 
defined as 

S aa (k,u) = £ |<gs|Sg(0)|n)| a % - (e n - E gs )), 

n 

where a denotes x,y,z, \n) are the excited states, and 
is the Fourier transform of Sf . As we work on fi- 
nite systems the SU (2) symmetry is not broken explicitly 



and S xx = S vv = S zz (in what follows the aa indices are 
discarded) . A straightforward calculation leads to the ex- 
pression 

S(K W) = Jjj^L/ |Uk+qWq - Wq«k+q| 2 <5(w - (w_ q + CJ k+q )), 

q 

(5) 

which satisfies the correct sum rule JJ2k a S aa (k,uj)duj = 
Ns(s + 1). 

As at the mean field level the triplet excitations are 
made of two spin-i free spinons a broad two spinon con- 
tinuum is expected. Nevertheless, as the 120° long range 
Neel order is developed there can be distinguished three 
distinct contributions in the spectrum. Following the in- 
terpretation of the spectra of [35] , it is instructive to split 
eq. ([5]) as 

by using the fact that u ± q = \v , q | ~ (^p-) 5 and w ± q ~ 
0. For k = ±Q, the spectrum is dominated by zero energy 
processes that create two spinons in the condensate. This 
gives rise to the magnetic Bragg peaks which, to leading 
order, behave as S±q u ~ Nm 2 S(uj). For k / ±Q, the 
spectrum is dominated by low energy processes that cre- 
ate one spinon in the condensate and another one in the 
normal fluid. This gives rise to a double peaked signal 
proportional to m, represented by 

ctsinq Td . |2 cv \ 

S K J = ^K + a-« k+ s I % -w k+ a) + 

m, l2r , > 

+ jl*« k _S+v k _fa| *(w-w k _|). 

Then, the shifted spinon dispersion w k± Q can be iden- 
tified with the low energy physical magnetic excitations. 
Finally, at high energy, the spectrum is dominated by the 
processes of creating two spinons in the normal fluid. This 
gives rise to a broad continuum represented by 

<Sk°™' -^yE'l Uk +1 t ' < l ~ U q U k+q| 2 <5(w ~ (W-q + Wk+q)), 

q 

where the prime means that sum goes over the triangular 
BZ except for q = ±^ or ±^ — k. 

In fig. [T]we have plotted eq. ([5} for the M point of the 
BZ (see inset of fig. [IJ. As noticed above, the low energy 
double peaked structure comes from S^™ 9 while the high 
energy tail corresponds to the continuum S^ ™*. In order 
to get the spectrum in the energy-momentum space we 
have plotted in fig. [2] the intensity curves of S(k,w) (eq. 
(O) along the path shown in the inset of fig. [1] The yel- 
low and red curves are the shifted spinon dispersion w, q 

of S^™ 9 while the blue zone corresponds to S^ ™'. In the 
figure we compare with the dispersion relations obtained 
with LSWT (solid line) and the recent series expansion 
calculations [T3] (points). At low energies the dispersion 
agrees quite well with LSWT and series expansions, be- 
ing the spectral weight mostly located around k ~ ±Q 
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Fig. 1: Dynamical structure factor, SQc, w), for momentum 
M — (§7r, ^tt). Inset: path of the triangular BZ along which 
the spectrum has been investigated. O=(0, 0), A=(ir,0), Q — 
(|tt,0), D = (27r,0),B = (7r,^=7r), and C = (§7T, -^=tt). w is 
measured in units of J. 



ABC O A Q D 

Fig. 2: Intensity curves for the dynamical structure factor, 
S(k,oj), calculated with the mean field Schwinger bosons the- 
ory within the two singlet scheme. Solid green line and blue 
points are the dispersion relations obtained with LSWT and 
series expansions [14], respectively. The path along the BZ is 
shown in the inset of fig. [T] 



(points Q and C). In this regime the physical excitations 
correspond to long range transverse distorsions of the lo- 
cal magnetization which are correctly described by both, 
LSWT and mean field Schwinger bosons. At higher ener- 
gies LSWT is not valid any more since the true spin ex- 
citations show a strong downward renormalization along 
with the appearance of roton-like minima (points). Re- 
markably, the mean field Schwinger boson theory predicts 
a non trivial redistribution of the spectral weight between 
the two spinon branches modulated by the form factor 
of eq. (0. The reconstructed dispersion, resulting from 
those pieces of spinon dispersion with the dominant spec- 
tral weight, reproduces quite well the series expansions re- 
sults. In particular, the crossing of the spinon dispersions 
at points B and D can be identified with the roton-like 
minima observed in series expansions. Regarding the in- 
terpretation of the roton minima, the singlet bond struc- 
ture of the Schwinger boson theory takes naturally into 
account the collinear spin fluctuations even in the pres- 
ence of the 120° Neel order of the THM. For instance, 
the roton minimum located at B can be interpreted as 
the development of magnetic correlations modulated by 
the magnetic wave vector (tt, "^ 7r ) which corresponds to 
certain collinear correlations pattern, while the other two 
non equivalent midpoints of the edges of the hexagonal BZ 
corresponds to different collinear fluctuations patterns. In 
fact, if these fluctuations are favoured by introducing spa- 
tially anisotropic or second neighbours exchange interac- 
tions the roton minima soften, giving rise to the new Gold- 
stone mode structure of the stabilized collinear ground 
state [2"Tp7]. 

Performing the frequency integration it is possible to 
analyze the relative weight of the two spinon continuum 
(blue zone of fig. [2]) to the static structure factor S(k) 
[2"8] . In fig. [3] we plot S(k) with diverging peaks located 
at the expected magnetic wave vectors ±Q (upper panel), 
along with the relative weight of the two spinon contin- 



uum, J / S(k)dLu (bottom panel). Interestingly, the 
contribution of the two spinon continuum to S(k) is ne- 
glegible around ±Q while outside their neighbourhood, 
and in particular at the roton position, the contribution 
to S(k) is about 40%. 

Density-density correlation functions. The small peak 
of fig. Q] leads to the remnant weak signal of fig. [5] 
which can be traced back to the local density fluctua- 
tion of Schwinger bosons. In fact, to describe the phys- 
ical Hilbcrt space of the spin operators the local con- 
straint of the Schwinger bosons must be satisfied exactly, 
Sf = + 1). Then, no fluctuations on the number of 

boson per site should be observed. However, since the con- 
straint is taken into account on average there are unphys- 
ical spin fluctuations in S(k, uj) coming from such density 
fluctuations. In order to identify them we have computed 
the density-density dynamical structure factor defined as 

N(k,u) = ^Kgs|n k (0)|n)| 2 ,5( W - (e n - E gs )), 

n 

where is the Fourier transform of the number of bosons 
per site, hi = J2 a "id 3 ™- A little of algebra leads to the 
expression 

N(k,Uj) = — ^ |uk+qVq + UqVk+q| 2 <S(w- (w_ q + W k+q )) , 

q 

(6) 

which is similar to eq. ([5]) , except to the plus sign within 
the form factor. If we split the two spinon contributions as 
N(k, u) = A£ 1 " 5 + iV£°"' it is easy to show that the main 
signal is located again at the shifted spinon dispersions 
U), o . But now, due to the different form factor, there is 

K + 2 

an important spectral weight transfer between such spinon 
dispersions. This is shown in fig. [H where we have plot- 
ted the intensity curves of N(k,u)) (eq. (J6))). It can be 
clearly observed that now the dominant signal is gapped 
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Fig. 3: Static structure factor (upper panel) and relative weight 
of the two spinon continuum, J /S(k)du, (bottom panel) 
along the same path of the BZ. 
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Fig. 4: Intensity curves for the density-density dynamical struc- 
ture factor, iV(k, uj), calculated within the mean field Schwinger 
bosons based on the two singlet scheme. The path along the 
BZ is shown in the inset of fig. [T] 



at Q and C points, while most of the spectral weight is 
located around k ~ 0. Such a soft mode can be iden- 
tified with a spurious tendency of the bosonic system to 
phase separation. Given the notable resemblance with the 
strong signal of N(k,uj), we suggest that the low energy 
weak signal of figs. [1] and [5] could be ascribed with the 
unphysical density fluctuation effects which we expect to 
disappear once they are projected out. For the unfrus- 
trated square lattice w^+(i i ) = ^k-(f ,f ) so both, the 
unphysical and the physical spin excitations, overlap in 
energy-momentum space, giving rise only to one low en- 
ergy band in 5(k,w) 0128]. 

Comparison with the one singlet scheme. So far we 
have found that the mean field Schwinger boson within 
the two singlet scheme reproduces quite well the series ex- 
pansions spectrum. It is also interesting to compare with 
the predictions of the one singlet scheme, since it is widely 
used in the literature. The first difference is the incorrect 
sum rule JJ2k a S aa (k, uj)duj = |iVs(s + I) which implies 
the well known | factor of Arovas and Auerbach [3] . Fur- 
thermore, in fig. [5] we have computed S(k, uj) after solving 
the corresponding self consistent equations for the param- 
eters Ag, and A. At very low energies the spectrum seems 
to be correct around points C, O and Q. However, at 
higher energies it is impossible to discern a reconstructed 
dispersion that fit the series expansion results along the 
whole path of the BZ, besides the factor about 3 in the 
energy scale. Therefore, we conclude that the two singlet 
scheme turns out the proper framework to describe cor- 
rectly the spectrum of the THM. Besides its quantitative 
accuracy, there are symmetry arguments that give further 
support to the two singlet scheme. In the literature, the 
one singlet scheme has been justified as the saddle point of 
a symplectic Sp(N) theory, originally adapted to extend 
previous large N works [3] to non bipartite lattices [5]. 
More recently, however, Flint and Coleman [55] demon- 
strated that if the Bij and Aij operators are kept the cor- 
responding large N extension preserves the time reversal 



properties of the spins, in contrast to the Sp(N) theory. 
Finally, it is worth to stress that the two singlet scheme is 
the basis of the Z 2 spin liquid theory, specially formulated 
to describe magnetically disordered phases [3"U] . 

Conclusions. — We have demonstrated that the sin- 
glet structure of the mean field ground state along with 
the fractional character of the spin excitations of the 
Schwinger boson theory take naturally into account the 
anomalous excitations of the spin-i triangular Heisenberg 
model recently observed |101fl4] . The appearance of the 
roton-like minima can be attributed to the tendency of the 
magnetic ground state to be correlated collinearly, even 
in the presence of 120° Neel order. By computing the 
density-density dynamical structure factor, and thanks to 
the series expansion results, we were able for the first time 
to discern, at the mean field level, between the physical 
and the spurious fluctuations coming from the relaxation 
of the local constraint. A further investigation within the 
context of the Schwinger boson theory reveals that the 
correct description of the spectrum depends crucially on 
the mean field decoupling. In particular, the two singlet 
scheme turns out more appropriate than the one singlet 
scheme. Based on the accurate description of the ground 
state static properties [11] (see table [1]) and in the light 
of the present results for the spectrum, we think that the 
bosonic spinon hypothesis should be considered seriously 
as an alternative viewpoint to interpret the physics of the 
triangular Heisenberg model. At the mean field level the 
triplet excitations consist of two spin-^ free spinons and, 
besides the low energy bands due to the onset of the long 
range order, there is a broad two spinon continuum, which 
could be related with the magnon decay found in the liter- 
ature [23]. In this sense, it would be important to improve 
the present mean field theory by deriving an effective inter- 
action between spinons resulting from 1 /N corrections or 
a better implementation of the constraint. We would ex- 
pect a picture of tightly bound spinons near the Goldstone 
modes while at high energies they would be weakly bound. 
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Fig. 5: Intensity curves for the dynamical structure factor cal- 
culated within the mean field Schwinger bosons based on the 
one singlet scheme. The path along the BZ is shown in the 
inset of fig. [1] Solid line and points are the same as in fig. [2] 

Work in this direction is in progress. Finally, we hope our 
present analysis in terms of bosonic spinons could help 
for a better understanding of the unconventional neutron 
scattering spectra of the CsiCuCl^ compound |17) . 

* * * 
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